clear
clc
close all
a = readmatrix('excel.xlsx','Range','A2:E18');
[m n] = size(a);
x0 = a(:,1:n-2);
y0 = a(:,n-1:n);
x0 = bbdesign(3);
r = corrcoef(x0);%计算相关系数矩阵
r = 1/4* (x0'*x0);
xd = zscore(x0);%对x0进行标准化处，均值化处理后，均值为0，标准差为1
yd = zscore(y0);%对y0进行标准化处理

[vec,lambda,rate] = pcacov(r);%计算相关系数矩阵的特征向量，特征值，各主成分的贡献率
contr = cumsum(rate);%累计求和
df = xd*vec;